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We report on both high-precision photometry from the MOST space 
telescope and ground-based spectroscopy of the triple system 5 Ori A con¬ 
sisting of a binary 09.5II+early-B (Aal and Aa2) with P = 5.7d, and a 
more distant tertiary (09 IV P > 400 yrs). This data was collected in 
concert with X-ray spectroscopy from the Chandra X-ray Observatory. 
Thanks to continuous coverage for 3 weeks, the MOST light curve reveals 
clear eclipses between Aal and Aa2 for the first time in non-phased data. 
From the spectroscopy we have a well constrained radial velocity curve 
of Aal. While we are unable to recover radial velocity variations of the 
secondary star, we are able to constrain several fundamental parameters 
of this system and determine an approximate mass of the primary using 
apsidal motion. We also detected second order modulations at 12 sepa¬ 
rate frequencies with spacings indicative of tidally influenced oscillations. 
These spacings have never been seen in a massive binary, making this 
system one of only a handful of such binaries which show evidence for 
tidally induced pulsations. 

Subject headings: binaries: close — binaries: eclipsing — stars: early- 
type — stars: individual (h Ori A) — stars: mass loss —stars: variables: 
general 


Introduction 


Massive st ars are rare in the Universe, due both to the low-mass favored IMF 
(|Salpeterl 1955) and their short life times. However, their high luminosities and 
mass-loss rates make them important to the ecology of the Universe. Underpinning 
this understanding is the observational determination of intrinsic stellar parameters 
for the most massive stars. However, this is not a simple proposition; except in rare 
cases where the distance is known reliably, it requires a detached, non-interacting 
binary system. We are aided by the fact that the binary fraction of O stars is nearly 
unity at birth, though thi s degrades so me as the stars evolve (e.g. mergers and 
dynamical interactions: iSana et al. 2014). Still, there are only m 50 O-star binaries 
in the Milky Way and Magellanic Clouds which meet the se cr iteria and for which 
individual stellar characteristics have been established (jGies 2012). This number may 
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increase in the future owing to advances in interferometry, enabling one to determine 
stellar m asses in systems s u ch as HD 150136 (03-3.5 V((f* ))+Q5.5-6 V((f))+06.5- 
7 V((f)), Sana et ah 2013 . Sanchez-Bermudez et. al 20131 ). Right now, the small 
number of systems spanning the O spectral type makes calibrating models difficult, 
but an excellent ef fort to de termine stellar parameters as a function of spectral type 
was made by Martins et ah ( 20051 ). 


Unfortunately, the problem of determining fundamental parameters becomes 
significantly worse for evolved O stars. It is extremely important to understand 
what effect leaving the main sequence has on the fundamental parameters of the 
star. However, there are currently only three such systems for which well-constrained 
masses have been derived, throug h the use of long baseline interferometry: C Ori A 


(09.2 Ibvar Nwk+09.5 Il-III(n), Sota et al. 2011, ISota et al.l 2014 . 


20131 ). HD 193322 (09.5Vnn+08.5III+B2 1 5y, 


Hummel et al. 


ten Brummelaar et al. 


2011) and MY 


Ser (07.5 III + 09.5 III-l-(O9.5-B0)IIII. ilbanoglu et al. 2013). While increasing the 
total number of massive star systems with well known parameters is worthwhile, 
it is difficult as the number of candidate O stars is quite small. A much more 
realistic goal is finding a non-interacting example system with minimal mass loss 
which can be used as a template for understanding the evolution of O stars. Mintaka 
(5 Ori), the westernmost belt star of Orion, by virtue of its proximity (d ~ 380 pc, 
Caballero fe Solano 2008 ) and brightness {V = 2.41, Morel fe Magnetat 1978 ). as 
well as encompassing the detached binary system Aa, is an excellent choice. 

The 5 Ori syste m has been known to contain multiple components for well over 
100 years flHa.rtma.nn et al. 1904). In fact, it is currently known to have at least 5 
stellar components (a full diagr ammed layout of this system and its relevant param¬ 
eters can be found in Fig. 1 of Harvin et ah ( 20021) ). The embedded triple system, 
especially the short period binary 5 Ori Aa, has received a significant amount of 
attention in the last 15 years. An a rtist’s renderi n g of the or bi ts and lay out of the 
6 Ori A can be found in Fig 1. of IShenar et al.l ((submitted]). lHarvin et al.l (12002 . 
hereafter H02), used spectra obtained with the International Ultraviolet Explorer 
to determine a putative secondary radial velocity curve. This, in combination with 
Hipparcos photometry, led to a controversial maximum mass for the O -star o f 11. 2 
Mq, less than half the expected value based on its spect ral type (see iMartins et al. 
2005). The system was revisited by (M a y er et al. 120101 . hereafter M10), in an at¬ 
tempt to understand this anomalous result. They worked under the assumption 
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that since the tertiary component contributed ~ 25% of the total light ( 

Perryman 

1997. 

Horch et al. 

2001, 

Mason et al. 

2009, 

Mafz Apcllaniz, J. 

2010, 

Tokovinin et al. 


2014), significantly more than the secondary, its presence might be able to explain 
the apparent discrepancy. They found that motion of the primary in combination 
with a wide tertiary component with spectral lines broader than that of the primary 
gave the false impression of orbital motion of the secondary. While M10 did improve 
the known orbital elements, the actual masses and radii are still poorly constrained. 
This is the result of a relatively weak secondary component as well as relatively poor 
light curve coverage and precision. 

This paper is part of a series which will attempt to characterize as rigorously 
as possible the triple system 5 Ori A, focusing on the_094311 primary Aal. This 
includes analysis of the X-ray sp ectrum bv ICorcoran et ahl (Isubmittedl . hereafter 
Paper I) and X-ray variabilit y bvlNichols et ahl (Isubmittedl . hereafter Paper II), as 
well as spectral modeling by Shenar et al. (Isubmittedl . hereafter Paper IV). In this 
paper we focus primarily on the optical variability. This includes an analysis of 
high-precision space-based photometry obtained with the MOST Space Telescope 
as well as simultaneous spectroscopic observations, described in detail in Section 2. 
In Section 3, we discuss the spectroscopic variability, including a derivation of the 
primary radial velocity (RV) curve. In Section 4, we use this RV curve in conjunction 
with the MOST light curve to provide three possible binary fits at varying values of 
the primary mass. In Section 5, we provide a fit to the primary mass using apsidal 
motion confirming our results from section 4. Finally, we discuss and characterize 
previously unknown photometric variability within this system and report frequency 
and period spacings consistent with non-linear interactions between tidally excited 
and stellar g-mode oscillations. While we are still are unable to disentangle the 
secondary spectrum and its kinematics, we provide a more complete model than has 
been achieved up to this point. 
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2. Observations 


2.1. High Precision Photometry from MOST 


Our optical photometry was obtained with the MOST microsatellite that houses 
a 15-cm Maksutov telescope through a custom broad-band filter covering 3500-7500 
A. The sun-synchronous polar orbit has a period of 101.4 minutes (/ = 14.20 d _1 ), 
which enables uninterrupted observations for up to eight weeks for targets in the con¬ 
tinuous viewing zone. A pre-launch summary of the mission is given bv !Walker et al 

fhoooh . 


PHASE 

0.0 0.5 1.0 1.5 2.0 2.5 3.0 3.5 4.0 



HJD-2450000 

Fig. 1.— MOST light curve (black) of 5 Ori taken in mid December 2012-early 
January 2013. Error bars are smaller than the size of the points. The red line is the 
combination of the binary fit (see Section 4) and that of the secondary variations 
(see Section 6). 




















5 Ori A was observed for roughly half of each MOST orbit from Dec. 2012 
through 7 Jan. 2013 with the Fabry-mode (see Fig. [TD at a cade nce of 40.4 s. These 
data were then extracted using the technique of IReegen et ah (120061 1. and show a 
point-to-point standard deviation of ~ 0.5 mniag. 


2.2. Optical Spectroscopy 

We initiated a professional-amateur campaign in order to obtain a large number 
of high-quality optical spectra simultaneous with our MOST and Ch andra campaigns. 
This resulted in more than 300 moderate resolution (R > 10000) spectra obtained 
over the 3 week period. The spectra were reduced by standard techniques utilizing 
bias, dark, and flat held images. Wavelength calibration was accomplished through 
comparison with emission-lamp spectra taken before and/or after the object spectra. 
Our main goal was originally to monitor the Ha line for wind variability and correlate 
any such emission with the X-ray spectra (Paper I, Paper II) while simultaneously 
obtaining a measure of the orbital motion from the He I 6678 transition. However, 
no significant Ha variability was seen, except for the radial velocity changes related 
to the orbital motion of Aal. Furthermore, Ha could be diluted by wind emission, 
making RV measurements questionable. Due to both the telluric absorption lines 
around Ha and the potential for wind emission, we focused on the He I 6678 transition 
for our radial velocity measurements. Details of the telescopes and spectrographs 
used are given in Table [H A typical S/N for any observation is m 100-150. 


3. The Spectroscopic Orbit 

The spectroscopic treatment of 5 Ori A from M10 and H02 are not in agreement. 
The M10 analysis assumed a tertiary contribution that was subtracted off, and left a 
single-lined binary, while the H02 analysis found evidence of the tertiary after iterat¬ 
ing to remove first the primary and then the secondary spectrum from the combined 
one. The resulting masses are highly discrepant, with M10 assuming a standard mass 
based on spectral type of M primary = 25M 0 , while H02 found M primary = 11.2 M & . 
Treatment of the tertiary may lead to severe errors in the inferred stellar parameters. 
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Paper IV presents a quantitative photospheric and wind model for each of 
the three stars in the 6 Ori A system. The results of iMason et al.l (120091 1 and 


Tokovinin et ah (2014) show that the tertiary orbits the binary Aa very slowly with a 


period of several hundred years, although the uncertainties in the tertiary’s orbit are 
very large. As the period is so long, we expect no measurable kinematic motion in 
the tertiary component during one month of observations. Therefore, if we take the 
model of the tertiary spectrum from the models of Paper IV, then we can subtract its 
contribution a priori from the spectra in order to best constrain the binary properties 
of Aal,2. This is similar to what was done by M10. If the tertiary is not taken into 
account, then our results are in agreement with H 02. However, the optical spectrum 
consists of ~ 25 % contribution from the tertiary flPerrvmanlll997l. iHorch et al.ll2001 . 


Mason et ah 2009), so even with our S/N ps 100-150, this contribution should always 


be seen in our spectra. 


We convolved the calculated spectrum of the tertiary to the resolution of our 
data, and can immediately see that it is the likely source of excess absorption. In 
Figure [21 we show two example profiles of the He I 6678 line (with large red and blue 
shifts from the primary’s orbital motion) along with the weighted contribution of the 
tertiary’s spectrum calculated by Paper IV. For the redshifted spectrum, we see that 
the red edge of the profile has an excess of absorption that matches the theoretical 
spectrum of the tertiary quite well. On the blue edge of the profiles, blending with 
He II 6683 becomes more problematic. While the tertiary seems to reproduce some 
of the absorption, most of blue excess comes from the primary star’s He II line. 

After we subtracted the weighted contribution from the tertiary star, we cross 
correlated our observations against the model of the primary calculated by Paper 
IV. From these velocities we calculated a single-lined orbit with the period derived 
by M10, as the length of our data sets was insufficient for determining a more pre¬ 
cise period. Most velocities had uncertainties on the order of ±2 km s^ 1 from this 
determination. Data for all the velocities obtained are given in Table 4. 

We examined the residuals after removing the primary star’s spectrum (both 
through shift-and-add and from model subtraction using the spectrum calculated by 
Paper IV). We found no residuals due to the secondary in the data at this point. 
This is consistent with the observational results of M10 as well as the results of the 
modeling of Paper IV. In particular, the modeling results predict that only 5.8% of 
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Fig. 2.— Two observed He I 6678 profiles (from HJD 2456279.99 showing the neg¬ 
ative velocity extremum, and 2456281.75 showing the positive velocity extremum). 
The He II 6683 absorption shown by a dashed red line represents a weighted non¬ 
moving spectrum of the tertiary calculated by Paper IV. 
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the light in this wavelength region comes from the secondary, so the noise associated 
with the weak features would dominate, as the individual spectra have a signal-to- 
noise on the order 100-150 on average. In order to reliably see the secondary with 
a S/N of 50, we suspect that the combined spectrum of the three component stars 
needs to have a S/N > 1000. 


4. A Full Binary Model 

A complete solution for a binary system requires a substantial amount of infor¬ 
mation. In general, it requires an eclipsing light curve and two radial velocity curves. 
Despite being studied numerous times, no one has been able to reliably determine 
all three of these quantities for 5 Ori Aa. This is primarily due to the faintness of 
the secondary with respect to both the primary and the tertiary. Despite the lack 
of a secondary RV curve, data have two significant advantages over the previous 
works of H02 and M10: high precision photometry along with simultaneous spec¬ 
troscopy. Since there is known apsidal motion in this system (see Section 5) having 
all photometry and spectra taken all within a month timespan mitigates this issue. 


4.1. Initial Setup and Modeling 


For simul ation and fitting of t he S Ori system we used the alpha release of 
PHOEBE 2.0 (IDegroote et al.ll2013l . iPrsa et al.ll20131) . an engine capable of simulat¬ 
ing stars, binaries, and concurrent ob servati onal data. Unlike the original PHOEBE, 
this is not built directly on the work (W ils on fe D ev inne v 19 71) and has been rewrit¬ 
ten to take into account many second-order effects. This is especially relevant when 
dealing with spaced-based data. 


Before this system was fit, some initial processing was required. As shown in Fig. 
[TJ the light curve of <5 Ori is composed of an eclipsing light curve with second order 
variations superimposed. These underlying variations cannot be easily or completely 
removed (see Section 6) so we must lessen their effects in order to get a more accurate 
representation of the binary variability. This can be done by first phasing on the 
binary period, which requires a precise ephemeris for the system. Since our data span 
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only 3 weeks (about 3.5 orbits), the period determined by M10 is much more precise 
than what we are able to obtain, so we adopt their value P = 5.732436 ± 0.000015 
d. For T 0 , min and T 0tPeriastron we adopted our own derived values. For T 0i p er iastron 
we determined a value from fits to the radial velocity curve of 2456295.674 ± 0.062 
ffJD. For To m i n we are actually able to obtain a value from O-C calculations of 
2456277.7897 ± 0.0235 ffJD. We should note that these numbers offer no specific 
improvement in the associated error over those derived in M10, and in the case of 
To.min the error is slightly worse. However, unlike M10 our method for determining 
To,min is independent of the binary fit. This is important as second order stellar 
effects (unknown in previous works) make determinations from the binary fit less 
reliable. In addition, apsidal motion is significant in this system, and we did not 
want it to affect our results. 

Using the revised ephemeris, the data were phased and then binned to provide a 
mild smoothing, which yields a template for the binary variations. In this template, 
the error is much smaller than the size of the points. However, incomplete removal 
of the secondary variations leads to scatter at the milli-mag level (see Fig. ED- 

In addition, modeling of stars and binary systems requires the presence of real¬ 
istic atmosphere models. Unfortunately, the Kurucz atmosphere tables are the only 
ones which have been incorporated into PHOEBE 2.0 at this stage and they are 
insufficient for O stars. Therefore, our best option is to use blackbody atmosphere 
models. In some cases this would not be adequate, but for O stars we sample the 
Raylcigh-Jeans tail at optical wavelengths, making the differences from a blackbody 
minimal. In our specific case, a comparison with non-LTE models calculated for the 
components of the system indicates a deviation of « 5% in the integrated flux across 
the MOST bandpass (Paper IV). Therefore, this assumption will not significantly 
affect our solution as all three stars have similar temperatures and deviations from 
a blackbody. 


4.2. Fitting and Solution 


PHOEBE 2.0 has several built-in fitting routines, but for the purposes of this 
project we decided to use the Monte Carlo Markov Chain method (MCMC). The spe- 
cihc implementation used by PHOEBE 2.0 is emcee (jForeman-Mackev et al.l 120131) . 
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PHASE (TO) 

Fig. 3.— The MOST light curve of 5 Ori phased on the orbital period. It is then 
binned to 0.0125 in phase to smooth out remaining secondary variations. The bottom 
x-axis is phased to time of minimum light and the top is phased to periastron. 
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Fig. 4.— The top two panels show the binned MOST light curve (left) and the 
radial velocity curve (right) overlaid with the low mass PHOEBE 2.0 fit (see Table 
[2]). In the bottom panels are the O-C residuals for the corresponding light and radial 
velocity curves. The intrinsic stellar variations that we are unable to completely 
remove are the main source of error. Only one model is plotted as the differences are 
nearly impossible to distinguish at this scale. 
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an affine invariant MCMC ensemble sampler based on Goodman & Wear e 020lOlh 


While this routine requires a good initial fit to converge to a solution in a reason¬ 
able number of iterations, it gives the most comprehensive analysis of both the fit 


of photometric time series (e.g. 


overview of this technique, see Gilks et ah ( 19961 ). 


Gregorv 

2005 

and 

Croll 

2006) 


However, this fitting method is not particularly useful without a full complement 
of constraints. Since we have no secondary RV curve, we must put an extra constraint 
on the mass. We know the spectral type of the primary star qui te we ll a nd ca n 
make a reasonable assumption that the mass is 25 ± 4M 0 flMartins et al.l 20051) . 
further supported by non-LTE modeling of the system by Paper IV. Unfortunately, 
as there are several viable sets of parameters, the fit was unable to converge to a 
single minimum. Therefore, we divided this mass range further into three separate 
ranges: low mass primary (21-24 M 0 ; LM) , medium mass primary (24-27 M 0 ; MM), 
and high mass primary (27-29 M 0 ; HM). These mass ranges were chosen based on 
previous results, as the system would converge to one of three possible solutions, 
with one falling in each bin. 


For each possible solution we iterated at least 800 times. While this may seem 
insufficient to converge to a true solution, we were not starting from a random point 
in the parameter space. Parameters were adjusted so that the overall shape and am¬ 
plitude were consistent with the observed light curve and RV curve before beginning 
MCMC. Indeed, after about 200 iterations the log probability function defining our 
goodness of fit is virtually constant. 

A representation of the quality of these fits for the LM model is shown in Fig. 
[4] with the upper panels showing the PHOEBE fit to the light curve (left) and RV 
curve (right). The lower panels show the corresponding residuals. While there are 
differences between the three models, they are impossible to distinguish on this scale 
and so only the LM model is shown. The only motivation for this model over the 
others is that its primary mass agrees best with what is derived from apsidal motion 
(see Section 5). One can see a comparison of the light curve residuals of each of these 
models in Fig. [5j This figure demonstrates the differences inherent in the residuals 
of each model. These differences are small but still significant, especially within 
the eclipses. However, it is impossible to determine quantitatively which solution is 
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best as the standard deviation of the residuals for each fit are the same to within 
« 2 x 10 -5 . It should also be noted that PHOEBE 2.0 is still undergoing rigorous 
testing, so while our models seem reasonable, we also simulated our solutions with the 
original PHOEBE (IPrsa & Z witter 200511 to ensure that the results were consistent. 


PHASE (PERIASTRON) 

- 0.2 0.0 0.2 0.4 0.6 0.8 1.0 



Fig. 5.— Shown is a comparison of the residuals of each fit: LM-MM (red), LM-HM 
(blue), and MM-HM (green). Each fit has a standard deviation of around 2.5 mmags. 
They are very nearly identical to each other in all places except in the eclipses. 

Table [3] shows the complete solutions for each model as well as comparisons to 
previous solutions by H02 and M10. While the spurious detection of the secondary 
RV biases the H02 results, there is good agreement between our work and that of 
M10. One point to note is the improvement made to the inclination angle. While 
each of the previous works had fits at substantially different inclinations from 67° - 
81°, our work shows that this inclination is quite firmly centered at 76° ±4° regardless 
of the fit. In addition all of the orbital parameters are well constrained. It is very 
difficult to say if the high, low, or medium mass fit is the best. The probability 
function is virtually identical for the three fits, and the only noticeable differences 
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are those for the masses and radii of the primary and secondary. In both the RV and 
light curve solutions the residuals are due more to the accuracy of the points than in 
any differences in the model (see Fig. [Hand Fig. 0). Therefore, we must allow for a 
family of solutions until we are able to retrieve a secondary RV curve. With no clear 
way to select a best fit model, we adopt the LM fit throughout this paper when we 
need to apply a model. 


5. Apsidal Motion 


In close binaries such as 6 Ori Aa, the deformations of the stars in the elliptical 
orbit cause the longitude of periastron (c o) to precess with time, which is referred 
to as apsidal motion. The amount of this precession depends o n the mass, inter- 
nal structure of the stars, and the orbital eccentricity. iBenvenuto et al.l (120021) and 
Ferrero et ahl (20131) made assumptions for structure constants of massive stars and 
then used the measurements of apsidal motion to calculate stellar masses. 

6 Ori A has been studied extensively for more than a century and shows signs 
of significant apsidal motion. The work of Harve y et ah (119871 1 compiled all deter¬ 
minations of the spectroscopic orbit. This is combined with more recent studies of 
H02 and M10 along with our own work to retrieve orbital information spanning over 
100 years. We use these derived orbital elements to calculate the apsidal motion, 
u, using a weighted linear least-squares fit to u as a function of time. The results 
shown in Figure [H] present a reasonable fit to the measurements and provide a value 
of uj — (1.45 ± 0.04)°yr _1 . 

With this measurement of the apsidal advance, we can make some assumptions 
about the internal structure constants of the primary star. This relics heavily on 
th e age of the syst em. The age of 5 Ori is likely best constrained by the study of 
Caballero fe Solanol (120081 ) who find 5 Ori Aal to be the brightest star in a small 
cluster. Their study involved deep imaging to understand low-mass stars, and their 
color-magnitude diagrams imply an age of ~ 5 Myr. In order to continue with the 
calculation we adopt our orbital parameters from §4.1. Combining these and the 
models of Benvenuto et al. (12002 ). we obtained the results shown in Fig. [3 which 
give the stellar mass as a function of age for our derived apsidal motion. With an 
age of 5 Myr, we see that the derived mass of the primary is ~ 22M 0 . This value 
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BESSELIAN YEAR 


Fig- 6.— Apsidal motion of <5 Ori. A weighted linear fit to the published values 
f Harvev et al. 1987 . H02, M10) yields a measurement of u — (1.45 ± 0.04)°yr _1 . 
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agrees with the spectral-type/mass (IMartins et ahl 20051) and our low mass model 
(Table [2]). 



Fig. 7.— Derived mass as a f un c tion of age from the apsidal motion and assumed 
spectroscopic parameters of 6 Ori. C aballero fe Solanol (120081 1 estimate the age of the 
cluster containing 5 Ori, to be 5 Myr, corresponding to a mass of « 22 M & consistent 
with our low mass model. The different lines represent the model determinations 
with the extremes of the apsidal motion from the errors (Fig. [6]). 


6. Photometric Time-Series Analysis 


It is clear from Fig. [T] that the dominant signal in the light curve comes from 
binarity. It is also apparent from the non-phase-locked variations that additional 
signals are present in 5 Ori A. The presence of such variation is not completely 
unexpected as iKholtvgin et ahl (120061 ) reported spectral variability in this system. 
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However, the variation timescales seen in FigJT] are much longer than the ~ 4h 
periodicities they report. 



Fig. 8.— The MOST light curve with the binary signal removed. The residuals 
show the presence of obvious secondary variations in the light curve. The black line 
shows the Period04 fit to the data. The lines show the location of the primary (red) 
and secondary (green) eclipses. 

Before addressing this variation, we first bin our light curve with the MOST 
orbital period (0.074 d). This allows us to mitigate the problem of light scattered 
into the MOST optics (mainly due to the bright Earth). As a consequence, we 
ignore variability on time scales shorter than a single orbit. Once this is done, we 
must remove the eclipsing binary signal. This variability has already been analyzed 
in the previous section, and hinders our ability to study the additional variations 
present. 

There are two options available for removal of the binary signal. We could use 
the binary template derived earlier (see Fig. |3]) or the LM model fit to this template. 
The decision of which binary light curve to use is not obvious, so the analysis was 
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done using both subtractions. The results are largely the same, and we chose the LM 
model fit subtraction as the amplitude of the significant peaks was higher and there 
were slightly more significant frequencies. This choice is supported by the agreement 
between the subsequent analysis and the full light curve (see Fig. [H red line). The 
residuals are shown in Fig. [8] (blue dots). 



Fig. 9.— Fourier transform of the binary subtracted residuals of 5 Ori A (black).The 
spectral window is plotted in green and shows no significant aliasing. The Fourier 
transform after pre-whitening is shown in blue. The red line at 1.49 millimag repre¬ 
sents the 4 a significance. The inset is a zoom in of the Fourier transform from 0-2 
c/d. Pre-whitened frequencies are labeled with arrows (see Section 6). Note that 
some frequencies are not immediately apparent without removal of larger amplitude 
peaks. 

A common method for finding and classifying periodic signals is known as pre¬ 
whitening. In this method peaks are detected in the Fourier transform and then 
a sinusoid with the corresponding period, amplitude, and phase is subtracted from 
the light curve. This is repeated, and a combination of sinusoids representing each 
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frequency are fitted and removed until there are no more significant peaks present. 
This method only works if the signals are sinusoidal (or mostly sinusoidal) and in 
the case of 5 Ori A, this is not obvious. 

The Fourier tran sform of t he 8 Ori A light curve is shown in Fig. [9] (black) 
using Period04 fjLenz fe Bregerl 2003). We apply pre-whitening to the transform 
and the end product is shown in blue. In order to determine which frequencies to 
inclu de as s ignificant, we use an empirically determined amplitude-to-noise ratio of 4 
flBreger et ah 1993 1 ). There are 12 frequencies (black arrows) above this significance 
threshold 4cr = 0.00149 (red line in Fig. [9]). The parameters associated with each 
frequency are given in Table [3J The fit of these frequencies to the residuals, which 
is_shown in Fig. [S] (black line), does a reasonable job of fitting the variations. While 
Kholtvgin et al.l (20061) attribute the variability they find to non-radial pulsations, 
we see no significant variability in the spectrum around 4 hrs (6 c/d). 

One possible explanation of extra periods is the presence of rotational moclula- 


Canticllo & Braithwaite 

(201 

1) lend some credibility to this idea. Moreover, the work 

of Ramiaramanantsoa et al. 

2014) 

Blomme et al. 

(2011 

), and 

Henrichs and Sudnik 

(2014 

) show evidence for the 

occurrence of this phenomenon among O stars. If rota- 


tional modulation is the root cause of the variation, then the largest frequencies would 
occur near the rotation period of the star. The vast majority of the remaining peaks 
could be explained as spurious signals due to t he variable nature of star spots ( e.g. 


lifetimes and number) as explored extensively by Ramiaramanantsqa et ah (2014) for 


the 07.5 III(n)((f)) star £ Per. In our case the period of the largest amplitude signal 
is 4.61 d. If we assume this this is the rotation period of <5 Ori Aal, then combined 
with the v sin i ps 130 kms -1 (Paper IV) and the assumption of alignment between 
the orbital and rotation axis, this yields a radius of ~ 12 Rq for the primary Aal. 
This is smaller than the value found in our models (see Tab. [4]) and closer to that of a 


luminosity class III star ([Martins et al.l 200511 as opposed to class II which 5 Ori Aal 


is known to be^ H oweve r, it is unclear what the uncertainty in radius would be in 


the models of Martins et ah 02 0 0 51 ) and so this peak remains roughly consistent with 


rotational modulation. Indeed, the phase folded light curves of the two strongest 
frequencies (see Fig. [TO]) show evidence of non-sinusoidal behavior. This could be a 
consequence of a highly non-sinusoidal signal such as rotational modulation. 
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Fig. 10.— Both plots show the phased residuals (small red points) overlaid with 
the binned (large black points) data to emphasise the shape. The top plot shows the 
largest amplitude peak (4.61 d). In the bottom plot this signal has been removed 
and the data is phased on the next highest frequency (2.49 d). 


While this hypothesis is plausible, similar variations can also be caused by pul¬ 
sations. While rotational modul ation would result in a forest of peaks that show no 
clear relationship to each other flRamiaramanantsoa et al.l 120141) . pulsations should 
show evidence of coherent spacings. All the periodicities we observe are several hours 
to days long, therefore th ey w ould fall in the g-mode regime for all three massive 
stars seen in this system ( Pamyatnvkh 19991 1. These modes have the unique quality 
that when they have the same l number they are equally spaced in period. Further 
explanation can be found in any asteroseismology review (e.g. Aerts et ajJ 2010. 


Now, we can check for period spacings by creating an echellc diagram where the pe¬ 
riod of each variation is phased modulo a given spacing and plotted. Equally spaced 
periods will lie roughly along a vertical line. The echelle diagram in Fig. fill shows 
that at least half of the frequencies share a common period spacing. The results of 
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0.05 0.10 0.15 0.20 0.25 

PERIOD MODULO 0.349 (D) 


0.30 0.35 


Fig. 11.— Echellc diagram phase-folded over a period of (0.349 d). Note the line of 
periods which fall around 0.05 indicating that they share a common period spacing. 


this analysis provide evidence for pulsations, but this evidence is no more convinc¬ 
ing than the idea of rotational modulation, ffowever, there also appear spacings in 
the frequency domain. Curiously, 10 of the 13 frequencies show a spacing at two 
or 3 times the orbital frequency (see Tab. SJ). This spacing suggests that there are 
ti dal interactions associa te d with these signals. This i s not unprecented as the work 
of Palate fc Rauw (2014), Koenigsberger et al. (2006) and Palate et al. ( 2013 1 also 
show some evidence of effects caused by tidal forces in massive binaries. 

This variability, though, is unlikely to arise from tidally exc it ed os cillati ons, 
which occur at integer ha rmonics of the orbital frequency ( Fuller et al.ll2012 . Burkart et al 
20121 . Wel sh et al. 1201 ll ). in contrast to the oscillations we observe. However, fre¬ 
quency spacing at multiples of the orbital frequency has been seen in several other 
close binary stars ( Hambleton et al~1l2013 . Borkovits et al. 2014 . Maceroni et al. 2014 1. 
and indicates that the pulsations are tidally “influenced.” Similar to the papers cited 
above we speculate that the frequency spacing at multiples of the orbital frequency 
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may arise from non-linear interactions between stellar oscillation modes and tidally 
excited oscillations. Moreover, we interpret the frequency spacing as evidence that 
the observed variability is generated by stellar oscillations rather than by magnetic 
activity or rotational modulation. 

The main caveat to this hypothesis is that our frequency resolution is limited to 
0.047 d~ x due to the short span (3 weeks) of the observations. This means that there 
are several spacings that could have been deemed significant. For instance F7 — Fll 
and F2 —Fll would both, within error, fit the same spacing because of how close the 
two peaks are in frequency. When this happens, we choose the peak that resulted 
in a spacing closest to a multiple of the orbital frequency. Additionally, the two 
frequencies which do not appear in Tab. GO F4 and F5, are actually split by the 
orbital frequency. However, since the uncertainty is half the actual spacing, we choose 
not to include it. Despite these issues, this frequency spacing is pervasive, appearing 
strongly between nearly every significant peak. While there is still uncertainty in 
this claim, it can be lifted with the addition of longer time baseline observations. 
Another problem is that even if these are tidally “influenced” oscillations, they still 
do not completely account for all the variation in this system. It is possible that 
these discrepancies could be due to amplitude modulation, though the ratio of period 
versus observation length makes this unlikely. It is also possible that there is also 
rotation modulation in this system from one of the components, just at a lower level 
than the pulsations. 

It is worth noting that the nature of the significant periodicities found in the 
lowest frequency range of the MOST data set is uncertain and could be a result 
of ’red noise’. Normally, this noise could be modeled and removed using Autore¬ 
gressive Moving Average techniques. However, our data is not stationary which is 
a requirement for this type of modeling. Attempts at making the data stationary 
disrupt the relevant signals; not only in our data, but also in the synthetic data using 
combinations of purely sinusoidal signals. Therefore, it is impossible to determine 
unequivocally the root cause of this variation without the addition of longer baseline 
photometric data (see Section 8). 


7. Discussion 


The system S Ori A remains an ideal candidate to define the intrinsic parameters 
of evolved 0 stars. The fundamental parameters of the primary, 5 Ori Aal, appear 
consistent with models from Paper IV as well as the results from modeling of the 
apsidal motion. Our analysis removes ambiguity in the derived stellar parameters 
and makes the primary more consistent with stellar models, without having to appeal 
to non-conservative mass loss or mass transfer that might play a role in the evolution 
of a close binary like 5 Ori Aa. 


Despite these consistencies, S Ori Aa has a rather larger eccentricity of « 0.114, 
which is unusual for a close O-star binary system with a moderate age of 5 Myrs. 
In such a situation, circularization of the or bit should already have occurred. We 
can check this using the procedure outlined in Zahn (119771 ) with one key assumption: 
knowledge of the stellar structure constant E 2 . This constant, which is dependent 
on the interior composition and structure of the star, requires extensive effort to 
calculate precisely, but can be approximated by the fractional core radius to the 
eighth power. We have no evolutionary models for 5 Ori Aa s o we have no way to 
calculate or estimate E 2 directly. We do know from Zahn et al. (1975J) that E 2 tends 
to increase as a function of mass, so we will take the value of E 2 « 3.49 x 10 -6 which 
is that of a 15 M 0 luminosity class V star. We choose this value as this is the largest 
mass star for which a value of E 2 was derived. This is clearly an underestimate of 
E 2 which is ameliorated slightly by the fact that 5 Ori Aa is a slightly expanded 
luminosity class II star (implying a smaller fractional core radius). It is also clear 
based on the inverse dependence of the circularization time with E 2 that the value we 
obtain will be a maximum. A quick calculation reveals a circularization timescale of 
« lMyr, which is well short of the age estimate that has been derived for this system. 
This result shows that our knowledge of this system is incomplete and could indicate 
that our assumed primary mass is incorrect. However, this discrepancy could also be 
caused by incorrect age calculations, strong interactions with the tertiary component, 
or issues with our knowledge of stellar structure in massive stars. At the same time, 
there is currently no evidence which strongly supports any of these hypotheses, so 
we are confident in the assumptions made. 


The most intriguing aspect of this system is the pronounced second order vari¬ 
ations in its light curve. Such variations have never been seen before in this system, 
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and the frequency spacings seem suggestive of tidally induced pulsations. If con¬ 
firmed, this would be an important discovery as these modes have been seen rarely 
if ever in massive stars. This could lead to asteroseismic modeling and increase our 
understanding of the interiors of both massive stars in binaries, and massive stars as 
a whol e, esp ecially considering the extremely high fraction of these stars in binaries 
f Sana et ah 2014 b 


While we have made significant progress towards understanding this system, it 
is clear that we are not at the level we need to be. Our ultimate goal is for § Ori to be 
a test for both evolutionary and structural models of massive stars. For this objective 
to be realized, we must determine fundamental parameters to high precision. For 
this reason it is imperative that we identify the presence of the secondary component 
within the spectrum. In addition, we must improve our frequency resolution within 
the Fourier spectrum. This will require a long time-series of high-precision pho¬ 
tometry. As a result, we could identify with confidence the source of the variation, 
which could possibly allow for asteroseismic modeling. A secondary consequence of 
this would be an improved binary light curve with significantly reduced scatter from 
which to model our system. This would also lead to improved values of fundamental 
parameters, and place constraints on the apsidal advance of the system. 


8. Future Work 

This analysis to understand 5 Ori, in conjunction with Papers I, II, and IV, has 
highlighted the need to obtain the spectral characteristics and RV curve of 5 Ori 
Aa2. We will use both the Hubble Space Telescope and Gemini-North to obtain high- 
resolution spectra of the binary at opposite quadratures while spatially separating 
the tertiary and observing its spectrum. These observations are scheduled for the 
upcoming observing season, and we also plan to obtain higher signal-to-noise opti¬ 
cal spectroscopy from the Observatoire de Mont Megantic. These observations will 
provide excellent constraints on the modeling of the system and reveal any unusual 
spectroscopic variability seen in the tertiary star. This will provide more constraints 
on the photometric variability seen with MOST . 

A better characterization of this variability as well as the binary parameters 
requires access to long time baseline photometry, which MOST is not able to pro- 





vide. The BRITE-Constellation project, however, is desi gned to provide 6 months 


of continuous coverage in both blue and red filters flWeiss et al. 


2014). BRITE- 


Constellation began initial data acquisition of the Orion constellation in October of 
2013. Since this was the initial data taken, the full 6 months of continuous coverage 
was not achieved. However, there are nearly three continuous months of data in both 
filters. These data should be released to the collaboration in early 2015. In addition, 
a second run with BRITE-Constellation began in Oct. 2014 that should provide a 
full 6-month time-series. 
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Table 1. 5 Ori Observing Log 


Observer 

N (Obs) 

Telescope 

Spectrograph 

CCD 

R 

Christian Buil 

4 

C9 (0.23 m) 

eshel 

ATIK460EX 

11,000 

Christian Buil 

6 

Cll (0.28 m) 

eshel 

ATIK460EX 

11,000 

Thierry Garrel 

70 

Meade LX200 14 (0.35 m) 

eshel 

SBIG ST10XME 

11,000 

Keith Graham 

2 

Meade LX200 12 (0.30 m) 

LHIRES III 

SBIG ST8Xme 

15,000 

Bernard Heathcote 

32 

Cll (0.28 m) 

LHIRES III 

ATIK314L+ 

13,000 

Thierry Lemoult 

5 

C14 (0.36 m) 

eshel 

ST8XME 

11,000 

Dong Li 

16 

Cll (0.28 m) 

LHIRES III 

QHYIMG2Pro 

15,000 

Benjamin Mauclaire 

8 

SCT 12 (0.30 m) 

LHIRES III 

KAF-1603ME 

15,000 

Mike Potter 

201 

C14 (0.36 m) 

LHIRES III 

SBIG ST8 

11,000 

Jose Ribeiro 

13 

C14 (0.36 m) 

LHIRES III 

SBIG ST10 

15,000 

Asiago 

4 

1.82 m 

REOSC echelle 


22,000 

Observatory UC Santa Martina 

33 

0.5 m 

PUCHEROS echelle 

FLI PL1001E 

20,000 

McDonald 

5 

2.7 m 

Cross-Dispersed 


60,000 

Nordic Optical Telescope 

4 

2.5 m 

FIES 


46,000 

Calar Alto 

3 

2.2 m 

CAFE Echelle 


65,000 


Table 2. Delta Orionis Fundamental Parameters Using Different Mass Constraints 


LM 


MM 


HM 


M10 


H02 


Mi(M 0 ) 

23.81 

24.20 

27.59 

25 

11.2T1.8 

M 2 (M 0 ) 

8.54 

8.55 

9.27 

9.91 

5.6±0.4 

til (Rq) 

15.12 

15.34 

16.08 

15.6 

13 

R2(Rq) 

5.00 

4.92 

5.17 

4.8 

5 

Tem(Kr 

30000 

30000 

30000 a 

30000 

33000 

T e f f2 (/f) 

24149i 3 7 4 21 

24039! 3 3 1 65 

r>QQQc628 

iOooo _297 

23835 


To,per+2450000(HJD) 

6295.674 ± 0.062 

b 6295.674 ± 0.062 b 

6295.674 ± 0.062 b 

4002.705 ±0.060 


I’o,min+2450000(HJD) 

6277.790 ± 0.024 

c 6277.790 ± 0.024° 

6277.790 ± 0.024° 

4001.966 ± 0.004 


P(days) 

5.732436“ 

5.732436“ 

5.732436“ 

5.732436 ± 0.000001 

5.723503 ± 0.000026 

UJ 

141.21191-^e 

141 43°- 75 ° 
_0. 7 5° 

141.0834 (_g 5 4 ° 80 

140±1.8° 


a(R & ) 

42.97_ 1 . 37 d 

43.14]_'g 5 i4 

44.86g.g° 15 

44±0.3 


e 

OH24^g° 0 9 1 2 01 

011 24.0-0092 
u.iiZ4_ 0 0101 

0.1130g.g o o 9 o 5 71 

0.0955±0.0069 

0.075±0.06 

i 

76.39 2 _T>°30 

77 23^'^° 

1 1 -4.10° 

76 ' 74 -i 5 16° 

73.6±0.3 ° 

77° 

'y(km s — ) 

15-51-I 3 76 

15.7lL? 6 i4 

15.34g.? 2 39 

21.7±0.5 


Ki (km s _1 ) 

96.02±0.60 

96.02±0.60 

96.02±0.60 

106.33T0.71 

94.9±0.6 


Note. — The errors quoted on the Low, Medium, and High models are equivalent to lcr errors, but are applicable to non-normal 
distributions as is often the case with MCMC. 

a value fixed during fitting, 
k value fit using only RV curve 

c value obtained from O-C calculations of primary minimum 

^Only one bound is given because the best fit is outside the lcr errors. 
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Table 3. Fitted Frequencies Due to 
Second-order Variations in 5 Ori A 


Name 

Frequency (c/d) 

Period (d) 

Amplitude (mag) 

Phase 

_F1 

0.402 ± 0.047 

2.49 ± 0.332 

0.0047 

± 0.0003 

0.289 ±0.011 

F2 

0.217 ±0.047 

4.614 ±1.284 

0.0038 

± 0.0004 

0.31 ± 0.016 

F3 

0.922 ± 0.047 

1.085 ± 0.059 

0.0034 

± 0.0003 

0.146 ±0.015 

FA 

0.155 ±0.047 

6.446 ±2.817 

0.0031 

± 0.0004 

0.525 ± 0.023 

F5 

0.331 ± 0.047 

3.023 ± 0.503 

0.0026 

± 0.0003 

0.066 ± 0.025 

F6 a 

0.034 ± 0.047 

29.221 ± 106.396 

0.0026 

± 0.0006 

0.421 ±0.041 

F7 

0.283 ± 0.047 

3.535 ±0.707 

0.0025 

± 0.0005 

0.985 ± 0.025 

FS 

0.99 ±0.047 

1.01 ±0.051 

0.0022 

± 0.0003 

0.855 ±0.021 

F9 

0.564 ± 0.047 

1.775 ±0.162 

0.0017 

± 0.0003 

0.06 ± 0.035 

F 10 

0.468 ± 0.047 

2.138 ±0.24 

0.0018 

± 0.0003 

0.63 ± 0.029 

Fll 

0.621 ± 0.047 

1.611 ±0.133 

0.0017 

± 0.0003 

0.59 ± 0.032 

F 12 

1.237 ± 0.047 

0.809 ± 0.032 

0.0018 

± 0.0003 

0.783 ± 0.025 

FI 3 

1.337 ± 0.047 

0.748 ±0.027 

0.0016 

± 0.0003 

0.985 ±0.03 

a This peak is likely an 

artifact due to a 

trend in 

the data. 

It is not 


considered real, but it is formally significant and included in the fit. 


Table 4. 


5 Ori A Orbital Frequency Spacings 


Frequencies Spacing (c/d) Significance 


F7-F11 

0.337 

2 forb 

F8-F13 

0.347 

2forb 

F2-F9 

0.347 

orb 

F9-F3 

0.358 

2 forb 

F8-F10 

0.522 

Sforb 

F1-F3 

0.520 

3 forb 



Table 4. Radial velocity curve 


HJD 

(d) 

RVel 

(krn/s) 

&RVel 

(km/s) 

6245.5063 

-30.0080 

10.00 

6261.4593 

-73.5121 

1.45 

6262.4548 

-49.8351 

1.41 

6263.4688 

35.0017 

5.00 

6264.4300 

98.0658 

1.30 


Note. — Table 4 is published in its entirety in the electronic edition of the As- 
trophysical Journal. A portion is shown here for guidance regarding its form and 
content. 







